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Abstract 



We study the Heisenberg antiferromagnetic chain with both dimerization 
and frustration. The classical ground state has three phases: a Neel phase, 
a spiral phase and a colinear phase. In each phase, we discuss a non-linear 
sigma model field theory governing the low energy excitations. We study the 
theory in the spiral phase in detail using the renormalization group. The 
field theory, based on an £0(3) matrix- valued field, becomes SO (3) x SO (3) 
and Lorentz invariant at long distances where the elementary excitation is 
analytically known to be a massive spin- 1/2 doublet. The field theory sup- 
ports Z 2 solitons which lead to a double degeneracy in the spectrum for 
half-integer spins (when there is no dimerization). 
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Antiferromagnets in low dimensions have been extensively studied in re- 
cent years, partly because of their possible relevance to high T c supercon- 
ductors and partly due to the variety of theoretical tools which have become 
available. The latter include non-linear sigma model (NLSM) field theories 
[1-6], Schwinger boson mean field theories 0, fermionic mean field theories 



series expansions ||, exact diagonalization of small systems [10], and 



the density matrix renormalization group (DMRG) method [11], [TI|. In one 
dimension, NLSM theories in particular have received special attention ever 
since Haldane [HJ conjectured that integer spin models would have a gap, 
contrary to the known solution for the spin-1/2 model, and this prediction 
was verified experimentally fl"3|| . 

In this Letter, we study a general Heisenberg spin chain with both dimer- 
ization (an alternation 5 of the nearest neighbor (nn) couplings) and frus- 
tration (a next-nearest neighbor (nnn) coupling J 2 ). Even classically (i.e., 
in the limit where the spin S — > oo), the system has a rich ground state 
'phase diagram', with three distinct phases, a Neel phase, a spiral phase and 
a colinear phase (defined below) ||14|| . For large but finite S, long wavelength 
fluctuations about the classical ground state can be described by non-linear 
field theories. These field theories are explicitly known in the Neel phase 
[H, and in the spiral phase (for 5 = 0) ^ ||. While the Neel phase has 
been extensively studied, various aspects like the ground state degeneracy 
and the low energy spectrum are not yet well understood in the spiral phase. 

We will first discuss the field theory in the Neel phase for arbitrary J 2 
and S. For the spiral phase, we show using a one-loop renormalization group 
(RG) analysis that the field theory flows to an 5*0(3) x SO(3) symmetric 



and Lorentz invariant theory with an analytically known spectrum |15| . We 



also discuss how the presence of Z 2 solitons (supported by the field theory) 
affects the ground state degeneracy and the low energy spectrum. Finally, 
we show that the field theory in the colinear phase is qualitatively similar to 
the one in the Neel phase. 

The Hamiltonian for the frustrated and dimerized spin chain is given by 

H = J x £ ( 1 + (-l) l 5 )S t ■ S m + J 2 £ S, ■ S l+2 , (1) 

i i 

where = S(S + l)h 2 , the coupling constants J\, J 2 > and the dimeriza- 
tion parameter S lies between and 1. Classically (for S — > oo), the ground 
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state is a coplanar configuration of the spins with energy per spin equal to 
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where 9\ is the angle between the spins S 2 i and S2M-1 and 6 2 is the angle 
between the spins S 2 i and S2i-i. Minimization of the classical energy with 
respect to the 9i yields the following phases. 

(i) Neel: This phase has 9\ = 9 2 = it and is stable for 1 — 5 2 > 4J 2 /Ji. 

(ii) Spiral: Here, the angles 9\ and 6 2 are given by 
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(3) 



where tt/2 < 9\ < 71 and < 9 2 < 71. This phase is stable for 1 — 5 2 < 
U 2 /J x < (1 - 5 2 )/5. 

(hi) Colinear: This phase (which needs both dimerization and frustration) is 
defined to have 9\=tx and 9 2 = 0. It is stable for (1 — 5 2 )/5 < 4J 2 / J\. 
These phases along with the phase boundaries are depicted in Fig. 1. 

We now study the spin wave spectrum about the ground state |TB[. A 
detailed analysis will be presented elsewhere [1?]. We only mention the qual- 
itative results here. In the Neel phase, we find two zero modes with equal 
velocities. In the spiral phase, we have three modes, two with the same 
velocity describing out-of-plane fluctuations and one with a higher velocity 
describing in-plane fluctuations. In the colinear phase, we get two zero modes 
with equal velocities just as in the Neel phase. The distinction between the 
three phases is also brought out in the behavior of the spin-spin correlation 
function S(q) in the classical limit. S(q) is peaked at q = ir in the Neel phase, 
at7r/2<g<7rin the spiral phase and at q = tt/2 in the colinear phase. 
Even for S = 1/2 and 1, DMRG studies have seen this feature of S(q) in 
the Neel and spiral phases [0]. The colinear region has not yet been probed 
numerically. 

The spin wave analysis is purely perturbative and is really not valid since 
there is no long-range order and no Goldstone modes in one dimension. To 
study non-perturbative aspects, we develop a NLSM to describe the low 
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energy modes. This is well-known in the Neel phase [fl], . The field variable 
is a unit vector <p and the Lagrangian density is given by 



Here c = 2J 1 aSVl — 5 2 — 4J 2 / ' J\ is the spin wave velocity (a is the lattice 



spacing) and g 2 = 2/(SJl — 5 2 — 4J 2 / ' J\) is the coupling constant. Note 
that large S corresponds to weak coupling. The third term in (4) is a topo- 
logical term with 9 = 27rS(l — 5). This field theory is gapless for 9 = n mod 
2tt with the correlation function falling off as a power at large separations, 
and is gapped otherwise. For the gapped theory, the correlations decay ex- 
ponentially with correlation length (, where ( is found from a one-loop RG 
calculation to be (/a = exp(2n/g 2 ). Hence ln(£/a) = nS^Jl — 5 2 — AJ 2 / J\. 
This is plotted in Fig. 2 for 5 = and 4J 2 / J\ < 1. 

Recently, the spiral phase has also been studied for 5 = [[5], The 
classical ground state has 9 1 = 9 2 = 9 = — J!/(4J 2 ). The field variable 
describing fluctuations about the classical ground state is an 50(3) matrix 
R(x,t) related to the spin variable at the i th site as (Sj) a = SJ2bRab n b > 
where a,b = 1,2,3 are the components along the x, y and z axis, and n is a 
unit vector given by 

x cos i9 + y sin i9 + ai . . 

I x cos i9 + y sin i9 + a£ | 

The unit vector rij describes the orientation of the i th spin in the classical 
ground state (assumed to lie in the x — y plane) and £ represents the local 
magnetization. The Hamiltonian in Eq. ([]]) can be expanded in terms of R 
and i and Taylor expanded upto second order in space-time derivatives to 
obtain a continuum field theory 0. The Lagrangian density is found to have 
an 5*0(3)^ x S0(2) R symmetry and can be parametrized as 

C = 1 tr(d t R T d t R P ) - °- tr(d x R T d x R P t ) , (6) 



where c = JiSa(l + AJ 2 / Ji)yl — J1/IQJ2, and P an d Pi are diagonal ma- 
trices with entries given by 

P = t J_ J_ 

\2gV 2gV g\ 2g 2 
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and P\ = 
The couplings are found to be 



1 
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and g\ = g 2 2 [ 1 + (1 - J X /2J 2 



(8) 



Perturbatively, there are three modes, one gapless mode with the velocity 
cg 2 / #4 and two gapless modes with the velocity cgi/g 3 . Note that the theory 
is not Lorentz invariant because gig^ ^ g2gs- However, the theory is sym- 
metric under SO(3)l x 5*0(2)^ where the S0(3)l rotations mix the rows of 
the matrix R and the S0(2)r rotations mix the first two columns. (To have 
the full S0(3)l x 50(3)^ symmetry, we need g\ = g 2 and g 3 = g 4 , i.e., both 
Pq and P\ proportional to the identity matrix.) The S0(3)l is the manifes- 
tation in the continuum theory of the spin symmetry of the original lattice 
model. The S0{2) R arises in the field theory because the ground state is 
planar, and the two out-of-plane modes are identical and can mix under an 
SO (2) rotation. The Lagrangian is also symmetric under the discrete sym- 
metry parity which transforms R(x) — > R(—x)P with P being the diagonal 
matrix (-1,1,-1). An important point to note is that there is no topological 
term present here (unlike the NLSM in the Neel phase) and hence, no ap- 
parent distinction between integer and half-integer spins. There is, however, 
a distinction due to solitons, as we will show later. 

At distances of the order of the lattice spacing a, the values of the cou- 
plings are given in Eq. (8). At larger distance scales /, the effective couplings 
giil) evolve according to the /^-functions (3{gi) = dgi/dy where y = ln(//a). 
We have computed the one-loop /^-functions using the background field for- 
malism ||18|| . (Note that since the theory is not Lorentz-invariant, geometric 
methods cannot be used to obtain the /3-functions |Q.) The /^-functions are 
given by 
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We have numerically investigated the flow of these couplings using the initial 
values gi(a) given in Eq. (8). We find that the couplings flow such that #1/(72 
and gz/g± approach 1, i.e., the theory flows towards SO(3)l x S0(3)r and 
Lorentz invariance. Finally, at some length scale (, the couplings blow up 
indicating that the system has become disordered. At one-loop, ( depends 
on J 2 / ' J\ but S can be scaled out. In Fig. 2, we show the numerical results 
for ln(£/a) versus J 2 / ' J\ for 4J 2 / J\ > 1- Note that as 4J 2 / ' J\ — * 1 from 
either side (the Neel phase for integer spin or the spiral phase for any spin), 
ln((/a) — > 0, i.e., the correlation length goes through a minimum. Since 
4J 2 /Ji = 1 separates the Neel and spiral phases, we may call it a disorder 
point. (For general 5, we have a disorder line 4J 2 /Ji + 5 2 = 1 and the 
correlation length is minimum on the disorder line separating the two gapped 
phases.) 

The spiral phase is therefore disordered for any spin 5* with a length scale 
(. Since the theory flows to the principal chiral model with 50(3)^ x S0(3)r 
invariance at long distances, we can read off its spectrum from the exact 
solution given in Ref. 
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The low energy spectrum consists of a massive 
doublet that transforms according to the spin-1/2 representation of SU(2). It 
would be interesting to verify this by numerical studies of the model. DMRG 



studies fil| , [L2[ of spin-1/2 and spin-1 chains have not seen these elementary 
excitations so far. It is likely that these excitations are created in pairs and a 
naive computation of the energy gap would only give the mass of a pair. To 
see them as individual excitations, it would be necessary to compute the wave 
function of an excited state and explicitly compute the local spin density as 
was done in Ref. to study a one magnon state in the Neel phase. 

Since the field theory is based on an S'0(3)-valued field R (x,t) and 
^1(5*0(3)) = Z 2 , it allows Z 2 solitons. The classical field configurations 
come in two distinct classes with soliton number equal to zero or one. If 
Ro(x,t) is a zero soliton configuration, then a one soliton configuration is 
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obtained as 

(cos$(x) sin#(x) 
-sm6(x) cosd(x) \ R (x,t) , (10) 
1 

where 6(x) goes from to 2n as x goes from — oo to +oo. (For convenience, 
we choose 6{x) = 2n — 9(—x), i.e., the twist is parity symmetric about the 
origin.) In terms of spins, this corresponds to progressively rotating the spins 
so that the spins at the right end of the chain are rotated by 2ir with respect to 
spins at the left end. Since the derivative d x 9 can be made vanishingly small, 
the difference in the energies of the configurations Rq(x, t) and R x {x, t) can 
be made arbitrarily small, and one might expect to see a double degeneracy 
in the spectrum. 

However, this classical continuum argument needs to be examined care- 
fully in the context of a quantum lattice model. Firstly, do Rq (x, t) and 
Ri(x,t) actually correspond to orthogonal quantum states? For the spin 
model, if the region of rotation is spread out over an odd number of sites, 
i.e., if the rotation operator is U — exp( 2r ^_ 1 5^=_ m (2n + 2m + 1)5^), then 
R J0 (x,t) and R\{x, t) have opposite parities because under parity, S* — > — Sf 
and U — > Uexp(i2n J2n=-m ^n)- Since the sum contains an odd number of 
spins, the term multiplying U is — 1 for half- integer spin and 1 for integer 
spin. Thus for half-integer spin, Rq(x, t) and Ri(x, t) are orthogonal and the 
argument for double degeneracy of the spectrum is valid. This is just a re- 



statement of the Lieb-Schultz-Mattis theorem [20|. For integer spin, i? (x, t) 
and R±(x, t) have the same parity and no conclusion can be drawn regarding 
the degeneracy of the spectrum. 

An alternative argument leading to a similar conclusion can be made fol- 
lowing Haldane . We consider a tunneling process between a zero soliton 
configuration R (x,t) and a one soliton configuration R, 1 (x,t). (We choose 
coplanar configurations for convenience). Such a tunneling process is not 
allowed in the continuum theory (which is why the solitons are topologically 
stable) because the configurations have to be smooth at all space-time points. 
But in the lattice theory, discontinuities at the level of the lattice spacing 
are allowed. In terms of spins, this tunneling can be brought about by turn- 
ing each spin S- ^ in configuration R (x,i) to the spin Sj in configuration 
R_i(x,t) by either a clockwise or an anticlockwise rotation. Assuming that 
the magnitude of the amplitude for the tunneling is the same (as we will show 
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below), the contribution of the two paths either add or cancel depending on 
whether the spin is integral or half-integral. This is easily seen through a 
Berry phase [22] calculation. The difference in the Berry phase of the two 
paths from Sjto is 2nS. Since the soliton involves an odd number of 
spins, the total Berry phase difference is mod 2tc if S is an integer and ir 
mod 2tt if S is half-integer. 

Now we have to check that the magnitudes of the amplitudes for tunneling 
are the same in both the cases. To see this, consider the pair of spins Sj 
and S^_l which need to be rotated to Sf^ and S^]. Since 9{x) = 2n — 9(—x), 
the magnitude of the amplitude for the clockwise rotation of to is 
matched by the magnitude of the amplitude for the anticlockwise rotation 
of SlI to S_]. Hence, for the pair of spins taken together, the magnitude of 
the amplitude for tunneling is the same for the clockwise and anticlockwise 
rotations. 

Thus, tunneling between soliton sectors is possible for integer S (thereby 
breaking the classical degeneracy and leading to a unique quantum ground 
state) but not for half-integer S (due to cancellations between pairs of paths). 
This agrees with the earlier Lieb-Schultz-Mattis argument. 

Although the NLSM model for the spiral phase was explicitly derived only 
for 5 = 0, we expect the same qualitative features to persist when 5^0, 
because the spin wave analysis shows that the classical ground state continues 
to be coplanar and there continue to be three zero modes (two with identical 
velocities and the third with a higher velocity [0). Hence we expect similar 
RG flows and a similar spectrum. However, the argument for the double 
degeneracy of the ground state for half-integer spins depends on parity being 
a good quantum number. When 5^0, parity no longer commutes with the 
Hamiltonian and the argument breaks down. This is in agreement with the 
DMRG studies Jl2| (for periodic chains) which show a unique ground state, 
both for integer and half-integer spins, for 5^0. For open chains, the ground 
state is sometimes degenerate due to end degrees of freedom. To incorporate 
such effects, one would have to study NLSM theories on open chains which 
is beyond the scope of this work. 

Finally, we examine small fluctuations in the colinear phase. The naive 
expectation is that the field theory would be an 0(3) NLSM, analogous to 
the Neel phase, since the classical ground state is colinear. We can show 



this explicitly for 5 = 1 which is called the Heisenberg ladder [2jJ. The 
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field theory in this limit can be derived using the classical periodicity under 
translation by four lattice sites, similar to the derivation in Ref. @ for the 
Neel phase. For two pairs of spins, we define 



(f)(x) = — , £(x 



2S y ' 2a 



and f(x) = S ^ 3 ~ s S ^ , H X ) = S ^+3 + S 4i+2 (n) 

We write the Hamiltonian in terms of the fields <fi and £, and Taylor ex- 
pand to second order in space-time derivatives to obtain the Lagrangian 



(§J) without a topological term. We now have c = AaSy J2 + Ji) and 

Q 2 = \ \J (J2 + Ji)/ -h- The absence of the topological term means that 
there is no difference between integer and half-integer spins and a gap exists 
in both cases. In fact, the NLSM predicts a gap for any finite inter-chain cou- 
pling, however small. This is in agreement with numerical work on coupled 
spin chains [E3~ . 



In conclusion, we emphasize that this is the first systematic field theoretic 
treatment of the general J\ — J2 — 5 model on a chain. It would be interesting 
to find an experimental system with sufficient frustration and dimerization to 
probe the colinear phase. This phase could also be studied using numerical 
techniques like DMRG. The field theoretic treatment of the spiral phase leads 
to the interesting possibility that the low energy excitations of integer spin 
models may be massive spin-1/2 objects. This again is a possibility which 
could be looked for experimentally or verified by numerical simulations. 
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Figure Captions 



1. Classical phase diagram of the J\ — J 2 — 5 spin chain. 

2. Plot of \n(C/a)/S versus J 2 /Ji for 5 = 0. For 4J 2 /Ji < 1, ln(C/a) is given 
by the one-loop RG of the 0(3) NLSM for integer spin. For AJ2/J1 > 1, 
ln(C/a) is given by the one-loop RG of the S0(3) L x S0(2) R NLSM. 
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